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Abstract. We show how spin-spin correlations, detected in a non-destructive way 
via spatially resolved quantum polarization spectroscopy, strongly characterize various 
phases realized in trapped ultracold fermionic atoms. Polarization degrees of freedom 
of the light couple to spatially resolved components of the atomic spin. In this way 
quantum fluctuations of matter are faithfully mapped onto those of light. In particular 
we demonstrate that quantum spin polarization spectroscopy provides a direct method 
to detect the Fulde-Ferrell-Larkin-Ovchinnikov phase realized in a one-dimensional 
imbalanced Fermi system. 



1. Introduction 

Condensation of fermionic pairs occurs in nature in a variety of systems, such as neutron 
stars, atomic nuclei, excitons in solids and superconducting materials. Ultracold Fermi 
gases have been in the last years at the forefront of research both theoretically and 
experimentally pQ. Indeed homonuclear mixtures of fermionic atoms in two different 
hyperfine (pseudo-spin) states offer the unprecedented advantage of a continuous 
tuning of the interspecies attractive interaction. Thanks to this unique feature 
recent experiments have spectacularly demonstrated high-temperature superfluidity of 
attractive fermions [2j [3j H], [5] for a large interval of values of the scattering length, 
spanning for the first time the crossover from Bardeen-Cooper-Schrieffer (BCS) pairing 
to the Bose-Einstein condensation (BEC) of bosonic molecules composed of two fermions 
(BCS-BEC crossover) [6]. The fate of the fermionic superfluid upon imbalancing the 
two spin species has been monitored in elongated traps, showing that conventional 
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pairing accompanied by segregation of the excess majority atoms persists over a large 
interval of imbalance values [7J |Sj. This finding stands in contrast to the expectation 
for one-dimensional systems, where exotic pairing with finite-momentum pairs of Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) type [9] is predicted to occur pH QU CGI H31 E] • In 
addition, ultracold heteronuclear fermionic mixtures [15] with controllable interactions 
have been very recently realized, as well Bose-Fermi mixtures in which the attractive 
interaction can create a Fermi gas of polar molecules [16J . 

One of the major difficulties confronted by recent experiments is the problem of 
faithful detection of correlations in the strongly correlated regime of the gas. For 
attractive fermionic gases, the best-developed probe for phase correlations to date is 
represented by the dynamical projection of fermionic pairs onto deeply bound molecules 
and the destructive time-of-flight measurement of the latter (2j [3, Ej. In the limit 
of deeply bound molecules, noise correlation analysis [I7J has been used as well to 
detect pairing correlations [18]. The other crucial probes used in experiments are 
phase-contrast imaging [7J E] which reveals only the local correlations between the two 
species; imaging of the vortex lattice induced by stirring, which reveals macroscopic 
phase coherence [5] ; and radio- frequency spectroscopy [31 [19] , which probes the binding 
energy of the fermion pairs. 

In this paper we focus on the fundamental insight that can be gained by shining 
polarized light onto the atoms and detecting the quantum fluctuations imprinted onto 
the light polarization by the atomic sample (quantum polarization spectroscopy, QPS) 
[20] . This non-destructive measurement gives direct access to spin-spin correlations in 
the atomic system [2~Tj . When the light shone on the atomic sample is a standing wave 
[22] , this type of measurement allows to directly probe the magnetic structure factor at 
the wavevector corresponding to the standing-wave period. Here we focus on the case of 
trapped attractive fermions, and we show that the measurement of the fluctuations of the 
light quadratures - which gives the information on spatially resolved spin correlations 
and fluctuations - is very sensitive to various aspects of the paired phases occurring in 
this system. In the case of a balanced Fermi gas undergoing a crossover from the BCS 
regime to the BEC regime, we show that the evolution of spin-spin fluctuations at a 
given wavevector exhibits directly the shrinking of the size of the pairs upon varying 
the scattering length. In the case of a one-dimensional imbalanced Fermi gas in an 
optical lattice, we show that FFLO pairing leaves an unambiguous fingerprint on spin- 
spin correlations. This signature enjoys the full robustness of the FFLO phase in ID 
[TOj [TTl [T2| [T3l [H] , and it persists also in presence of a parabolic trapping potential. 

The structure of the paper is as follows. In section [2] we introduce the atom-light 
interfaces and the interaction Hamiltonian and show how the atomic spin fluctuations 
are mapped onto the fluctuations of light. Section [3] deals with balanced Fermi mixtures, 
and the fate of spin-spin correlations across the BEC-BCS crossover. Section Hlfocusses 
on an imbalanced Fermi mixture in a one-dimensional optical lattice, and demonstrates 
the FFLO fingerprint on spin-spin correlations. Conclusions are discussed in Section |5j 
Before proceeding further, let us stress that the method we propose here to unveal spin- 
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Figure 1. (a) Setup for a single beam with intensity J(r) propagating in direction 8; 
(b) Two plane waves counterpropagating in direction 9, giving rise to a standing wave 
configuration in the propagation direction. The second one is obtained by a mirror 
reflection of the first one. 



spin correlations in ultracold atomic gases is not restricted to atomic samples. It could 
be extended, e.g., to study the spin physics of fermionic molecules created in Bose- Fermi 
mixtures. Fermionic molecules might posses very large electric dipole moments which 
modifies substantially the properties of the strongly interacting regime with respect to 
the case of short-range interactions. QPS could hence help to diagnose the effects of 
dipole interactions on the spin correlations in such a regime. 



2. Atom-light interfaces 

A thorough derivation of the atom-light interafaces that result from the propagation of 
polarized light in an atomic sample can be found in [23] . The effective dipole interaction 
between an atom with spin J and a linearly polarized off-resonant light propagating in 
an arbitrary direction reads 

H int = -a J dr a i ph l at + s z J z - -J= (a ] z a + + a}_a z j J+ - -J= {a\h- + a\a^j J_ 

Here, a^(r,t) and a(r,t) denote the creation and annihilation electric field operators, 
s z (r,t) = | ^a+a + — b)_d^ is the Stokes operator, and l a t( p h) denotes the total density 
of atoms (photons). The constant coupling a = a cj\ 2 h/(AS7i) where 7 is the excited 
state linewidth, A the wavelength of the probing laser, A the cross section of the probing 
laser overlapping with the atomic sample, 5 is the detuning and ao is the standard AC 
Stark shift (for quantum number F = 3/2, ao = 3). For a two-component Fermi system 
the atomic spin J can be defined in terms of the field operators (r, t) in the usual 
manner 

Jz = \($h-i>\j>i), J x = l(J + + J-), J y = i(i + -J_), (2) 

where J + (r,t) = tptyi, J_(r,t) = t/wj. 

We consider the different probing configurations shown in Fig. [TJ They lead to an 
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effective interaction Hamiltonian of the form^ 

= -a [ drA(r)s z J P (3) 



where Jp = J z cos 9 — J y sin 9 and the intensity A(r) = I(r) for a tilted beam shown in 
fig. [I] (a). Probing setup (b) leads to A(r) = 2 cos 2 [krp + <f)], where rp = z cos 9 — y sin 9 
and (p is the spatial phase of the probing standing wave. 

In general both the atomic spin vectors and the Stokes operators are functions of 
space and time coordinates. Neglecting retardation effects, we can integrate over the 
duration of the probe pulse (of the order of [is) and define a macroscopic Stokes operator 
S a = J dts a where a — x,y, z. We consider light initially polarized in the x- direction 
(S x ) = Np/2 and (S y ) = (S z ) = 0, where Np is the photon number. Heisenberg 
equations of motion for the atomic spin lead to J a (r,t) = J a (r) \/a to first order in 
time, while the spin in the direction of the probe Jp is conserved to all orders. Due 
to the atom-light interaction, the Stokes operator performs a Faraday rotation in the 
plane perpendicular to propagation 

d r S y = -aJ P {r)^A{r). (4) 

Integrating equation (J4J, we find that the light quadrature, defined by the Stokes 
operators S y , X s = ^/2/NpS y reads after propagation through the sample 

(^D = (XT) --£=( dr{J P (r))A(r), (5) 



where k 2 = o 2 NpNa/2, being Na the total number of atoms. The fluctuations read 

(8Xr) 2 = \ + ^J drA ( r ) J dr f A(r')(5Jp(r)SJp(r')), (6) 

where 5Jp = Jp — (Jp). Thus fluctuations of the light quadrature (polarization 
fluctuations) after crossing the atomic sample contain the photon shot noise plus 
a contribution proportional to second order correlations of the atomic spins. 
Experimentally, spin-spin correlations can be detected if their contribution in the 
polarization fluctuations is larger that the photon shot noise, which for a coherent initial 
source corresponds to (5X™) 2 = |. To this end, the relevant parameter is k 2 = rja 
where a is the resonant optical depth of the sample and rj is the spontaneous emission 
probability [20J, [22] . It can be shown that the optimal signal is obtained when r] is tuned 
to rj op t = l/y/2a and k 2 = yct/2 [30J. BEC clouds have typical optical depths of the 
order of a few hundreds, for which one obtains quantum fluctuations imprinted on light 
by the spin fluctuations which are significantly bigger than the photon shot noise. 

Technically it might be challenging to fix the spatial phase of the probing light 
with respect to the trapping potential of the atoms in standing wave probing 

| Here and in the following equations, the coordinate labels (x,y,z) of the Stokes operators are 
associated with a reference frame in which the beam propagates along the z direction. 
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configuration (b). We assume that the average of the signal over various shots of the 
experiment leads to averaging over 0. Under phase averaging, Eq. ([6]) reduces to 

{8Xf\k)f = i + [4S m (0) + S rn (2k) + S m (-2k)} . (7) 
Here we have introduced the magnetic structure factor 

S m {k) = \\ dr \ dr'e^ r '\5J P (r) 5J P (r')) (8) 

where V is the volume of the system and k is the momentum of the probing standing 
wave. If the Hamiltonian of the system conserves the total spin along the probing 
direction P, as it will be the case in the following examples, then S m (0) = 0. Moreover, 
due to the commutation of 5Jp(r at different locations, S(— 2k) = S(2k). 

In atoms are trapped in a one dimensional optical lattice as considered in Section 
IU the atomic spin field (e.g. along the z direction) is localized at lattice sites, and it is 
conveniently expressed as 

J z (r) = ^^2\w(r -ri)\ 2 rhi (9) 

i 

where £V runs over the sites of the optical lattice and rhi = {h^ — n.jj)/2 where 
is the occupation of fermions with spin a at lattice site i. w(r — r») is the Wannier 
function at lattice site i; in the following we will approximate it with a 5-function for 
simplicity. If the probing standing wave is not phase-locked spatially to the one creating 
the optical lattice, averaging over the relative phase between the two leads to Eq. (J7J 
which contains the lattice magnetic structure factor: 

S m (k) = M^l^Hn-r 3 ) _ . (10 ) 

id 

where w is the Fourier transform of the Wannier function, and L is the linear dimension 
of the (i-dimensional lattice. In the following we consider for simplicity w(r) ~ 8(r), 
and consequently we neglect the fc-dependence of w. 

Hence the method proposed here gives direct access to the magnetic structure 
factor of the atomic sample. This piece of information is crucial in detecting the onset 
of spin-spin correlations in the strongly interacting phases of trapped atomic samples. 
The (pseudo)spin degree of freedom of spin-S* atoms is encoded in 2$ + 1 internal 
hyperfine states which are populated in the atomic sample. A fundamental example 
of spin correlations is the antiferromagnetic phase emerging at low temperatures in 
a Mott insulator of spin- (1/2) bosons or fermions in an optical lattice [2H [25]. More 
exotic spin states appear in spin-1 bosons, including dimerized states in one dimensional 
optical lattices [26]. Moreover, recent proposals envision the dipolar coupling of a 
S = 1/2 pseudospin degree of freedom of molecules dressed with microwaves [27]: the 
spin couplings realized in this way give rise to novel magnetic phases including, e.g., 
topological quantum order. 
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3. Three dimensional balanced Fermionic superfiuids 

We begin by considering superfluidity in an homogeneous two-component balanced 
3D Fermi system with two-body contact attractive interactions. The presence of 
interactions leads to pairing between the different spins described by the pairing order 
parameter A(r) = \g\{^(r)iliUr)) , where the interaction strength g oc 2/cpa s , a s is the 
two-body scattering length, and kp is the Fermi vector. For an homogeneous balanced 
3D gas (J(r)) = and the spin fluctuations are the same along all directions. Using 
Wick's theorem, the fluctuations in J z read 

A(5J z (r)5J z (r')) = £«M P M P " O - K&(r)Vfe(r'))| a ) - 2|<^(r)^(r')>| 2 . (11) 

<?=n 

The BCS formalism [6] describes the BCS-BEC crossover, characterized by different 
values of the scattering length l/k F a s , in terms of the Bogoliubov amplitudes u \X v \r) — 

5 ( 1 ± ] ^ ~~ ^= ) where the chemical potential a and superfluid gap A are obtained 

from the simultaneous solution of the gap and number equations [28] : 

-^ = (/i 2 + A 2 ) 1 /4p ( X ) (12) 
k F a 

- = Kli 2 + A 2 ) 1/4 Pi /2 (x) + (fx 2 + A 2 ) 3 / 4 P_ 1/2 (x) (13) 

where x = — fx/ (fx 2 + A 2 ) 1 / 2 , A = A/ep, fx = fx/ep and P v (x) are the Legendre functions 
of the first kind. 

We consider two types of probing configurations: a running Gaussian beam 
[Fig. Q](a)] and a phase-averaged standing-wave configuration [Fig. Hfb)]. A Gaussian 
probing beam has an amplitude profile A(r) = e~\ ri -\ 2 l a2 , where r±_ are the coordinates 
perpendicular to the propagation direction. When the light is macroscopically polarized 
in the x-direction propagating through the sample, it yields (X s out ) = (Xf 1 ) and 

1 ^ L -^/^e-^ 2 (|/ TT (r)| 2 + |/ iT (r)| 2 ) 



(**r(k))' 



(14) 



where 



J TT( r ) = TTTZyi I dk k2 v l M kr ) 



o 



In( r ) = j^y J dk k' 2 v k u* k jo(kr), (15) 

and j is the Bessel function. 

For a homogeneous probing set-up the total spin fluctuations are zero [29], while a 
finite Gaussian probe shows a finite value of the spin fluctuations. 

Results in Fig. [2] show how the polarization fluctuations decrease with decreasing 
pair size along the BCS-BEC crossover. These results are in agreement with those of 
Ref. [31], and they reveal that the quantum fluctuations of the atomic spins imprinted 
on the light polarization are significantly suppressed when the characteristic pair size 
becomes smaller than the Gaussian beam waist a. This is due to the simple fact that 
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Figure 2. Spin fluctuations ((SX° ut ) 2 - in the BCS (a < 0)-BEC crossover 

for a Gaussian probing beam. Sizes kpo — 0.5 (solid) , kpo = 1 (dashed), kpo = 5 
(dotted) and k F a = 10 (dot-dashed). Here V = L 3 with Lk F = 100. 





(b) 



Figure 3. Probing schemes of spin fluctuations in paired fermionic gases, (a) Gaussian 
laser beam, (b) Standing wave. 



an s-wave pair is in a total spin singlet, and hence it does not contribute to the spin 
fluctuations of the atoms illuminated by the laser beam when it is fully contained within 
the beam waist. A sizable quantum noise imprinted in the polarization corresponds to 
the optimal situation in which the pair size is larger than or comparable to the beam 
waist, as sketched in Fig. [3](a). Yet a fundamental remark is necessary: in order to 
obtain a total excursion of order 0(1) in the noise signal along the crossover, one needs 
to focus the laser to the experimentally challenging waists a ~ 1 — 10 kp 1 . In fact, for 
typical sample densities of order n ~ 10 13 atoms/cm 3 , one has kp 1 = (3tt 

2 n )-V3 _ 150 

nm. 
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Figure 4. Spin fluctuations ((5X° ut ) 2 — for different momenta of a phase- 

averaged probing standing wave k in the BCS (a < 0)-BEC crossover. This is also 
proportional to the magnetic structure factor S m (2k). Different lines correspond to 
k/kp — 0.25, 0.5, 1, 2, 3 in ascending order. 

On the other hand, probing the system with a standing wave, as proposed in 
the previous section, gives high spatial resolution without the need of focusing a 
laser over prohibitively small length scales. Fig. H] shows the magnetic structure 
factor S m (2k)V /Na, Eq. (JSj) , at fixed wavevector k, obtained via the phase-averaged 
polarization fluctuations of two counter-propagating beams at wavevector k, as a 
function of the product (/c^a s ) _1 . Strongly correlated spin fluctuations, recorded by 
the magnetic structure factor, are present over length scales related to the pair size, 
while the inter-pair spin correlations are vanishing. Hence the magnetic structure factor 
is very strongly affected by the shrinking of the pairs controlled by the scattering length, 
which gradually pushes the correlation length of spin fluctuations below the probe 
wavelength 2n/k (Fig. E](b)). The light wavevectors k/kp = 0.25, 0.5, 1, 2, and 3, 
considered in Fig. HI correspond to wavelengths ~ 3770, 1885, 942, 471, and 314 nm 
respectively. The strongest excursion on the magnetic structure factor along the BCS- 
BEC crossover is observed for the largest wavevectors, which correspond to conventional 
laser wavelengths. The other wavelengths can be obtained effectively by crossing the 
counter-propagating beams at an angle, as we will discuss in a forthcoming publication. 

The inflection point of the magnetic structure factor S m (2k) as a function of (kpci) -1 
contains the information about the characteristic length ^ beyond which the spatial 
average of spin fluctuations imprinted on the light polarization is strongly suppressed. 
The various scans of S m (2k) at fixed k for varying {kpa) -1 allow to extract the (kpa)^ 1 
location of the inflection point, and to reconstruct how the associated wavelength 
£i = 2ir/k depends upon (fcpa) -1 . This dependence is shown in Fig. |5j It is remarkable 
to see that the length £j follows the same behavior of the Pippard coherence length £ p 
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Figure 5. Solid line shows the Pippard coherent length (£ p = 1/ An) in log-scale while 
the dashed line shows the corresponding length = 27r/fc where fco is the inflection 
point of the magnetic structure factor S m (2k). 

(up to a scaling factor) on the BCS side of the crossover, while it starts decreasing much 
faster around the unitarity limit. In this respect, it would be desirable to go beyond 
mean-field theory in the description of the strongly interacting regime, and to extend 
to that regime the comparison of the spatial features contained in the structure factor 
with the characteristic pair size. 

The above analysis can be applied to BCS pairs as well as to molecules. Therefore 
QPS based on standing waves appears as a promising method to probe the internal 
structure of s-wave Feshbach molecules. Further analysis of the QPS signal in the case 
of, e.g., p-wave or <i-wave pairing and Feshbach molecules will be the subject of future 
investigations. 

4. Imbalanced Fermions in one dimension. Fingerprint of the FFLO phase 

The possibility of exotic pairing in imbalanced fermionic mixtures has been recently 
the subject of a very intense research. In particular the mismatch between the Fermi 
momenta kp^ and kpi of the two spin species can lead to the appearance of FFLO pairs 
with finite momentum Q = \kp^ — kpi\. The experiments carried out so far seem to 
rule out this possibility in three dimensions, where phase separation of the gas into a 
balanced mixture (with conventional pairing) and a fully polarized gas of the remaining 
majority atoms is observed [7J[8]. On the other hand, analytical and numerical studies 
have rigorously proved the occurrence of a highly stable FFLO phase in one dimension 

[mi eh na nan]. 

In the following we concentrate on the case of one-dimensional imbalanced attractive 
fermions in an optical lattice, described by the ID attractive Hubbard model in a 
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parabolic trap: 

H = -t {4,A+i,a + h.c.J -U^2 n iT n a (16) 

2,(7 1 

where c\ a creates a fermion with spin a =f| at site i , n ia = c\ a Ci^, t is the hopping 
matrix and U is the on-site attractive interaction. 

We investigate this model by means of numerically exact quantum Monte Carlo 
simulations based on the Stochastic Series Expansion algorithm [32] in the canonical 
ensemble [33]. The temperature is set to reproduce the T = properties of the system. 
As pointed out in Section [2] the method of quantum polarization spectroscopy gives 
access to the magnetic structure factor Eq. ffTQ]) . which is the quantity we calculate 
numerically with quantum Monte Carlo. 

It is well known that, for a bipartite lattice, the attractive Hubbard model can 
be mapped onto the repulsive one via a particle-hole transformation on one of the two 
species, cL = (—1)* cVf, ct| = c 'li- The fillings of the d fermions are n'^ — 1 — rif and 
n\ = Hi, and accordingly the Fermi wavevectors are transformed as k' F + = ir — kp-^ and 
k' F < = kFi- m particular the spin-spin correlation function of the attractive model maps 
onto the density- density correlation function of the repulsive one 

C mm {r) = (rhim i+r ) - (rhi) (m i+r ) = - ((^-n- +r ) - n' 2 ) (17) 

where n! { = hL + h^t . The density-density correlation function for the repulsive Hubbard 
model has been investigated extensively in the past via Bethe Ansatz, and its general 
form is of the type 

^ — , cos( aik'v f r + 3ik'r?:r) 

CUr) = m + r) -n' 2 = J2 B i r A* (18) 

i 

where the various terms appearing in the sum are imposed by selection rules [34). In 
particular it is found that (ai,j3i) = (2,0), (0,2), (2,2), .... To the best of our knowledge, 
the amplitudes Ai of the various contributions have not been determined for the case of 
a spin-imbalanced Hubbard model at an arbitrary filling. 

The spin-spin correlation function for the attractive imbalanced ID Hubbard model 
is shown in Fig. [6] for various values of the attraction and of the imbalance. Our 
numerical findings are all consistent with the expression 

CW0 = \ c L( r) « M -L + , 2 ros(2 ^' r) = M -L + , 2 sm&H, (19) 

which, in the symmetrized form C mm {r) + C mm (L — r), provides excellent fits to the 
finite-size numerical data. 

This shows that the spin-spin correlation function of the attractive model is directly 
sensitive to the Fermi momentum mismatch Q at which FFLO pairing occurs. A sketchy 
picture on the relationship between spin-spin correlations and FFLO pairing in the 
strong pairing limit is offered in Fig. [71 The FFLO state is effectively a mixture of 
bound pairs and of unbound majority atoms, which are mutually hardcore repulsive 
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Figure 6. Spin-spin correlations in the ID attractive Hubbard model with population 
imbalance for different polarization values and interaction strengths. The simulation 
data refer to a system size of L = 80. To improve readability, in both panels some sets 
are multiplied by a constant (indicated in the set legends). The dashed lines are fits 
to the symmetrized form of Eq. (|19p . All results show consistency with the estimates 
Ai « 1.8 and A2 w 2 — 2.2, although the precision on the fitting coefficients does not 
allow us to extract the systematic dependence of the exponents on the polarization 
and interaction parameters. 



due to Pauli exclusion principle. Infinite repulsion in ID leads to algebraically decaying 
charge- density- wave (CDW) correlations, namely the system displays strong fluctuations 
towards a local CDW state which very roughly corresponds to an equally spaced 
arrangement of pairs and excess | particles [35] . From the point of view of the spin 
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Figure 7. Sketch of the short-range properties of an imbalanced mixture of strongly 
attractive ID fermions. a) The two-component gas displays strong fluctuations towards 
a short-range crystalline-like arrangement with the majority spins at a maximum 
mutual distance l/fif, and similarly for the minority spins, b) The strongly attractive 
interaction leads to pairing fluctuations, which set the two short-range crystalline 
arrangement in phase spatially. From the magnetic point of view, pairing effectively 
erases some of the itinerant spins in the system, so that the resulting magnetic texture 
at short range has a characteristic length « (n-f —n^)~ 1 and a characteristic associated 
wavevector 2tt{tl^ — n-f) = 2Q. 



texture, considering bound pairs as spinless objects leaves out an algebraically decaying 
spin-density-wave arrangement which has a characteristic wavevector 27r(rif — n\) = 2Q. 
Hence this argument shows that the 2Q-modulation of spin-spin correlations is a direct 
fingerprint of pairing in an imbalanced mixture. 

The 2Q modulation of spin-spin correlations translates into a pronounced kink at 
k = 2Q in the magnetic structure factor, as shown in Fig. [8] for various values of the 
attraction and of the imbalance. This kink is shown to be the only relevant feature in 
S m (k), and to be most pronounced for small imbalance. Fig. [8] also shows a comparison 
with the case of two non-interacting spin species U = 0, which is exactly solvable. In 
that case the magnetic structure factor reads (for k F ^ > kp\) 



\k\ 

S m (k) = 

7T 



\k\ + 2k Fl 
2^ 

kp-\ + k 



Fl 



71 



for \k\ < 2kpi 

for 2kp^ < \k\ < 2kp^ 

for 2k < \k\ < 7i . 



(20) 



This expression clearly exhibits two independent kinks at 2k p-^ and 2k As shown 
in Fig. [8], upon increasing the attraction \U\ among the two species, the two kinks of 
the non-interacting case disappear gradually, while the kink at 2Q appears clearly for 
\U\/t > 4. The occurrence of a single kink at 2Q in the interacting model is hence a 
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Figure 8. Magnetic structure factor in the ID attractive Hubbard model with 
population imbalance for different polarization values and interaction strengths. The 
arrows in the lower panel mark the location 2Q/tt. 



direct consequence of pairing. 

To make full contact with experiments, we also add a parabolic trap to the 
Hamiltonian Eq. ( 1T61) 

W = H + V t Y,^-L/2) 2 n i>ai (21) 

i.CT 

where L is the size of the lattice. We use similar parameters to those of Ref. [12J, 
which shows clear evidence of FFLO pairing in terms of off-diagonal correlators. Fig. [9] 
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Figure 9. Upper panel. Magnetic structure factor in the trapped ID attractive 
Hubbard model with population imbalance for different polarization values and 
interaction strengths. Middle panel. Momentum derivative of the magnetic structure 
factor. Lower panel. Corresponding density imbalance profiles, An = n-\ — n^. The 
dashed lines mark the correspondence between the central density imbalance in the 
trap, the position of the kink in the structure factor, and the corresponding jump in 
the structure factor derivative. 



remarkably shows that the kink feature in the magnetic structure factor survives to the 
presence of a trap. In particular, the presence of a kink is even better evidenced by the 
momentum derivative dS m (k)/dk, which correspondingly exhibits a marked jump. The 
persistence of this feature in a trap is intimately connected with the fundamental fact 
that the density imbalance An = —n^ is almost constant over a significant portion of 
the trap center (although the density profiles of the two species independently are not 
as flat). Given that the location of the kink is only sensitive to the imbalance Q = it An, 
the persistence of the same An over a large portion of the cloud protects the kink from 
smearing. Conversely, the position of the kink can be regarded as an efficient measure of 
the density imbalance in the trap center, which is not in principle known a priori. This 
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is a valuable alternative to the direct in-situ phase contrast imaging recently applied to 
polarized Fermi gases [36] for the measurement of the density profiles of both species. 

5. Conclusions 

We have shown how spin-spin correlations of an attractive Fermi gas can be detected 
in a non-destructive way using spatially resolved quantum polarization spectroscopy 
(QPS). The atomic spins couple to the polarization degree of freedom of light and 
they imprint their correlations on the quantum fluctuations of the polarization, which 
can be measured using homodyne detection. When a standing wave is shone on the 
atomic sample, the measured signal allows a high resolution of spin-spin correlations in 
momentum space at the wavector of the standing wave. 

Spin-spin correlations are shown to strongly characterize the superfluid phases of 
fermionic systems with attractive interactions. In a three-dimensional spin-balanced 
system the spatial structure of spin-spin correlations depends strongly on the pair size 
and hence it evolves strongly along the crossover from the BCS to the BEC regime: as 
a consequence, QPS is able to reveal the evolution of the characteristic pair size as a 
function of the scattering length. In a one-dimensional system with spin imbalance, the 
magnetic structure factor, recorded by the QPS signal, shows a kink at the difference 
between the Fermi vectors of the two spin species, providing a direct signature of 
finite-momentum Fulde-Ferrell-Larkin-Ovchinnikov pairing. The proposed experimental 
technique is most promising to detect the spin structure of exotic pairs and molecules, 
including e.g. p-wave and <i-wave pairing, and to detect the magnetic phases which 
can be potentially realized by strongly correlated atoms and molecules loaded in optical 
lattices. 
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